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Abstract 



We present a mathematical framework for generating thick domain wall 
solutions to the coupled Einstein-scalar field equations which are (locally) 
plane symmetric. This approach leads naturally to two broad classes of wall- 
like solutions. The two classes include all previously known thick domain 
walls. Although one of these classes is static and the other dynamic, the cor- 
responding Einstein-scalar equations share the same mathematical structure 
independent of the assumption of any reflection symmetry. We also exhibit a 
class of thick static domain wall spacetimes with different asymptotic vacua. 
Our analyses of particle motion in such spacetimes raises the interesting pos- 
sibility that static domain walls will possess a unique experimental signature. 
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I. INTRODUCTION 



Spontaneous symmetry breaking in the early universe can produce domain walls (some- 
times called cosmic membranes). In the simplest model, a domain wall is formed when a 
discrete symmetry is broken with the resulting discrete set of vacua in different regions of 
space-time 0. The study of domain walls has a long history dating back to the work of 
Zel'dovich et al 0, ( see also || for a recent review by Cvetic and Soleng ). Domain walls 
produced at high temperatures are cosmologically problematic since they produce large-scale 
anisotropies which violate current experimental bounds However, domain walls that 

are formed after the decoupling of the microwave background do not produce large-scale 
anisotropies @. These late-time domain walls are produced at low temperatures ( ~ 10 -2 
eV) and therefore have characteristic widths on the order of 1 MPc and low-energy densities, 
thus the name soft domain walls. Soft domain walls have large density gradients and thus 
are a potential source of the density fluctuations that are necessary for the formation of 
the large scale structure seen in the universe. It is also interesting to note that the energy- 
momentum tensor for domain walls violates the strong energy condition which results in a 
gravitational field that is, on average, repulsive |||9|]. 

In this paper, we consider thick domain walls that arise from a single real scalar field with 
self-interaction. The associated Einstein-Scalar field equations admit accelerated spherically 
symmetric (and hence, locally plane symmetric) wall-like solutions. Unlike the previous 
investigations §||l(| of other thick domain walls, we do not assume any reflection symmetry 
of the corresponding spacetime, and our analysis of the Einstein equations leads to two 
broad classes of wall-like solutions which include previously known domain walls. Although 
non-reflection symmetric domain walls have been previously investigated by Cvetic et al and 
by Jensen and Soleng P,pT]-p!5|] those walls were thin domain walls. In this paper our interest 
is in thick domain walls. Furthermore, we find that relaxing the assumption of reflection 
symmetry leads to a family of static plane symmetric domain walls. This result is consistent 
with the Dolgov-Khripovich theorem ||16|. The Dolgov-Khripovich theorem rules out static 
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domain walls provided the vacuum energy density is positive and the metric has the same 
asymptotic limit on both sides of the wall. By violating the metric condition we are able to 
produce explicit solutions for static domain walls. Linet JHJ has shown that one can also 
produce a static wall by using a negative vacuum energy density. 

These walls have the unusual property that their asymptotic vacua differ from each other 
intrinsically. 

In section II, we introduce a simple model of domain walls, and discuss the general 
properties of the associated energy-momentum tensor with particular emphasis on various 
energy conditions. Without reflection symmetry we obtain, in section III, a new class of 
domain walls which are (locally) plane symmetric. In section IV, the existence of static 
domain wall solutions is established. Finally, in section V, we present a detailed analysis 
of the nontrivial particle motion in domain wall spacetimes. In the Appendix, we prove 
two lemmas, which are essential for understanding the properties of static and accelerated 
domain walls. 

II. STRUCTURE OF THE ENERGY MOMENTUM TENSOR 

Following H we consider a real self -interacting scalar field on a space-time (M 4 , g). The 
action for the scalar field is given by 

A[0] = -/ A/ {iA 2 d0(d0)+mM0)}^l (2-1) 

where A and m set the energy scales, * is the Hodge dual operator and the vector field d(f> 
is the metric dual to d<p. [Throughout this paper, we work with units where h = c = ks = 1 
and G = (1.2 x 10 19 GeV)~ 2 . In these units <p and the 'potential' V{4>) are dimensionless.] 
The stress energy tensor due to <fi is given by 

T = A 2 d0 <g> d<p - |^A 2 d0 (df) + m 4 V(0) j g. (2.2) 

We assume that the potential V(<f)) is positive. This ensures that the energy density 

T(Z, Z) = ^ | [Zm 2 + E [^(0)] 2 } + ™V(0) (2.3) 
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is positive. Here {X , Xi, X 2 , X 3 } form a local orthonormal basis and Z = X is any observer 
with g(Z, Z) = — 1. 

Since we are looking for domain wall solutions we require that V(<f>) has two degenerate 
global minimum at — /i± with V(p±) = 0. The boundary between these regions is 
the domain wall where the potential V((f>) > and has a maximum at = /x • The 
field is a smooth function on spacetime which connects the vacua /i±. Thus, the (ideal) 
boundary must be a timelike hypersurface in (M 4 ,g), defined by = /io. For such field 
configurations the gradient field, d0, must be normal to the boundary hypersurface, and 
hence it is spacelike. The stress energy tensor can now be written as 

T = -p (g - N ® N) + v (N ® N) (2.4) 

— r / — — m V 2 

where the unit normal is defined by iV = d0/ g ( d0, d0J , and the functions p and v 
are given by 

p = lA 2 d0 (d$) + mV(0) (2.5) 

v = ]-\ 2 d(j) (d$) - mV(0). (2.6) 

Since V(4>) > and d0 (d0^ > 0, we have 

p > (2.7) 
p + u>0 (2.8) 

In order to search for a domain wall configuration of the scalar field as a solution to 
Einstein-Scalar equations 

we have made two physically reasonable assumptions: 

(i) \/(0) > 

(ii) is spacelike, which means that g ab d a (f)d b (f) > 0. 
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These two assumptions lead to the following properties of the stress tensor T: 

1. J a b satisfies the weak-energy condition: 

For an observer, u a (where u a u a = —1), we have 

u a u a T ab = X 2 (u a d a (f>)(u b d b( f>) + { Y9 ab d a <f>d b <f> + m A V{<t>)} > 

2. Tab satisfies the dominant-energy condition: 

Defining v a := T ab u , we have v a = \ 2 (d a 4>)(d b 4>)ub — pu a , which implies 

VaV a = -p 2 - 2m 4 AV(0)( Ma 9» 2 < 

3. T a b violates the strong-energy condition, u a u a J a b > — |^ r (T), for every causal vector u a : 
To see this we choose an observer u a with u a u a = — 1 and u a d a (f) = 0. Then, we have 

n a u a lab = {^g ab d a (t>d b( t> + m 4 V(0)} 

~tr(J) = {y^ a 0^0 + 2m 4 V(0)} 

Hence, u a u a T a b < — \tr(T), violating the strong-energy condition for the chosen causal 
vector field u a . 

4. Einstein's equations imply the equation of motion for the scalar field 0: 

= d a <PT ab ; b = (409 a 0){A 2 V 6 VV - ™ 4 ^^} 

Here, (<9 a 0<9 a 0) ^ since (<9 a 0) is spacelike, and hence Einstein's equations lead to the 
scalar-field equation 

A 2 V 6 VV-- 4 ^ = 

d(j> 

III. ACCELERATED SPHERICALLY SYMMETRIC DOMAIN WALLS 

In order to determine the spacetime due to the stress-energy tensor (2.2) we must solve 
the Einstein equations 

Ric-iR£ = 87rGT (3.1) 



where Ric and R are the Ricci tensor and scalar curvature respectively. To look for a 
spherically symmetric solution we choose a chart (t, x, 9, ip) and assume that the metric 
admits three Killing vector fields 

Kx = sin ipdg + cot 9 cos ipd^ 
K 2 = — cos ipd e + cot 9 sin ipd^ 

K 3 = (3.2) 
We consider a spherically symmetric metric of the following form: 

g = f(t, x) {-dt <g> dt + dx <g> drr} + h 2 (t, x) [d9 <g> d9 + sin 2 6> ® d^} 

(3.3) 

where the functions f(t,x) and h(t,x) are not assumed to be spatially reflection symmetric. 
Now, taking Lie-derivatives of Einstein's equations, (3.1), with respect to a Killing vector 
field K, we find 

L K T = 

U(traceT) = (3.4) 

which, together with the equations (2.4)-(2.6), gives K(p) = 0, K(z/) = 0, and hence K(0) = 
0. Then from (2.2), we have dg(f> = d^cp = 0. Since d<p is spacelike, is assumed to be 
t-independent in the chart (t,x,9,ip). Hence d(p = (d x <p)dx, and the spacelike unit vector 
field, N = d(f)/[g(d(j)id(f))]^ , normal to the membrane is given by N = (l/f)d x . Then, the 
membrane stress tensor (2.4) can be written as 

T = -p(-e° <g> e° + e 2 ® e 2 + e 3 ® e 3 ) + ve l ® e 1 (3.5) 

where we have defined the coframes by 

e° = fdt ; e 1 = fdx ; e 2 = hd6 ; e 3 = hsm9dip (3.6) 

and the corresponding orthonormal frames are 

X = (l/f)d t ; X, = (l/f)d x ; X 2 = (l/h)d 6 ; X 3 = (l//isin0)c^ (3.7) 
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For complete specification of the metric [(3.3)] inside a spacetime region dominated by a 
thick (cosmic) membrane, it is necessary to solve for the functions {/, h} from the Einstein 
equations 

P a = 87iG{J ah -(T/2)g ab }e b (3.8) 

where T is given by (3.5), and T = trace[T] = (—3p + u). Now, computing the Ricci 1-forms 
P a = Ric(X a , Xb)e b with respect to the orthonormal basis [(3.6)-(3.7)], we find 

p = Ae° + Be 1 ; P x = Be + Ce 1 ; P 2 = De 2 ; P 3 = De 3 (3.9) 

where {A, B,C, D} are functions of {t,x} [Henceforth, a dot and a prime denote partial 
differentiations with respect to t and x\. Now, from (3.5), (3.8)-(3.9), it is seen that B = 
R\c(X , Xi) = 0, and hence, in terms of the coeffecients {/, h}, we have 



A simple solution to (3.10) can be given by 

f = f(x) ; h = c(t)f(x) (3.11) 
Setting B = [(3.10)] in (3.9) and using (3.11) in the Einstein equations [(3.8)], we find 

fA ee d x ( f j) + 2 (L ] 2 - 2 Q = - 47 rG(p - z,)/ 2 



n,. 2(f 

/2 C =_ 3 ^^] =4 7rG(3p + z/)/ 2 



-fD ee 9 X f^l + 2 (T| 2 - f ^ - f^V - f-V = -4^( P - z,)/ 2 (3.12) 



// v/y Vc/ Vc/ Vc/ 

From the first and the third equations in (3.12), we have (c/c) = (1/c 2 ) + (c/c) 2 which has 
a solution 

c(i) = icosh(fct) (3.13) 
where A; > 0. Using (3.13) and defining u = f'/f, the Einstein equations (3.12) reduce to 

2u' + u 2 - k 2 = -{MGf 2 )p (3.14) 



- 3(w 2 - k 2 ) = -(MGf 2 )v (3.15) 

As a consequence of (3.14)-(3.15) or V • T = 0, the equation of motion for the scalar field <fi 
is given by 

It is clear from (3.5) that p and v represent energy density and pressure (normal to the 
wall) with respect to the orthonormal basis in (3.6)-(3.7) : J(X ,X ) = p; T(Xi,Xi) = v, 
and by (3.14)-(3.15), p = p{x) and v = v(x) in the chart {t,x,y,z}. [Here we remark that 
components of the stress tensor [(3.5)] are not, in general, time-independent with respect 
to the coordinate basis. For example, J(d y ,d y ) = — pexp(kt) by (3. 5), (3. 11) and (3.13)]. 
Furthermore, v — implies (f'/f) 2 = k 2 , and hence p = - leading to an empty spacetime 
(vacuum). Thus the pressure, (v), inside the membranes can not be zero. We now show 
that the nowhere zero smooth function v(x) < based on the constraint (p + v) > 0, and 
additional physical assumptions that / is bounded and nowhere zero, and the stress tensor 
[(3.5)] tends to zero in the limit \x\ — > oo giving rise to asymptotic vacua. 
3.1 Negative 'Pressure' (is) : 

From equations (3.12) and (3.14)-(3.16), we have 

3u' = -ATrG(3p + u)f 2 (3.17) 
v' = -3(p + v)u (3.18) 

From the equation (3.17), it follows that u'{x) < 0, and hence, u has at most one zero. 
Also, for the nowhere zero smooth function, u, there exists an x such that u'(x) = (See 
Appendix, Lemma 2, for a proof of this statement). Then, (3.18) implies that each of the 
functions u and v 1 has a unique zero at x — x. Furthermore, differentiating (3.18) and using 
(3.17), we have v"[x) = — (p + v)u'(x) > - which means that the unique extremum of v 
over R is a minimum at x — x. However, v satisfies the asymptotic condition lim^i^oo v — 0, 
and hence v must be negative. As a consequence of v < 0, we have p > \u\ since p > and 
(p + v) > 0. 



3.2 The Equation of State V = v(p)' : 

In order to solve the Einstein-scalar equations (3. 14)- (3. 15), we need an equation of state 
v = v(p) subject to the constraints v < 0, p > \v\, lim^^^ v = and lim^^^p = 0. 
First, we observe that u'(x) < \fx, and hence, from (3.15), lim x ^± OQ u(x) = =fk. From 
these properties of u(x), it follows that the non-zero smooth functions {u 2 — k 2 ) and u' are 
both negative; each of them has a unique minimum and tends to zero in the limit \x\ — > oo. 
Using (3.14) and (3.15), we also have 

. 2 U ' {X ] V 2 = -«*) (3-19) 

where the nowhere zero smooth function e(x) is given by 

Sp(x) 



> 1 (3.20) 



\U{X)\ 

since p(x) > \v(x)\. In fact, the equations (3.19) and (3.20) lead to an equation of state 
v = — VqP where (3 — Vq) /2vq = e(x) > 1, and hence, v% < 1. It is important to note that 
e(x) is dimensionless according to our choice of units [specified at the beginning of section 
II]. Thus, e(x) is, in fact, a function of the dimensionless argument, kx, where the dimension 
of k is given by inverse of length. 

To determine the classes of admissible e(x) leading to domain wall configurations of 4>(x), 
we recall from section II [and the discussion below (3.3)] that lim x ^± OQ <p(x) = p±, where 
the constants p + and represent two degenerate vacua with V(p±) = 0. However, from 
the equations (2.5)-(2.6), (3.3) and (3.14)-(3.15), and the existence of asymptotic limits of 
0, we have 

(j)'(x) 2 = (2/X 2 8nG) [e(x) - 1] [k 2 - u(xf 
lim <f>'{x) = (3.21) 

|x|— »oo 

where lim^i^oo [k 2 — u(x) 2 } = 0. Thus, any e(x) from the class of bounded smooth functions 
such that 1 < e(x) < e , will satisfy the condition (3.21). 

One can also choose e(x) to be unbounded. In this case, we have another family of ad- 
missible smooth functions satisfying (3.21) provided, asymptotically, [k 2 — u(x) 2 } approches 



zero faster than e(x) tends to infinity. To see this, we integrate (3.19) to obtain u(x): 

u{x) = — /ctanh (^k j e(x) dx + c^j 

where c is an integration constant. Inserting the above expression for u(x) in [k 2 — u(x) 2 ] , 
we find 

'k 2 - u(x) 2 

The first equation in (3.21) now gives 

<P\x) 2 = {2k 2 /X 2 87iG) [e(x) - 1] sech 2 (k J e(x) dx + c) 

Since sech 2 {k J e(x) dx + c) tends to zero exponentially as \x\ — > oo, we may take e(x) to 
be, for example, a polynomial such as 2 + 3k 2 x 2 for which (3.20) and (3.21) are satisfied. 

This suggests a procedure for solving the Einstein equations, (3.14)-(3.15), with a choice 
of e(x) from these two large classes of functions leading to domain walls. In the following 
subsection, we will discuss solutions due to bounded e(x). 
3.3 Solutions to the Einstein Equations : 

A simple solution to the Einstein equations, (3.14) and (3.15), is obtained by choosing 
e(x) = e > 1 - which trivially belongs to the class of smooth bounded functions. Hence, 
such a choice ensures, through (3.21), the existence of a domain wall solution to the Einstein 
equations. Integrating (3.19) and using u = /'//, we find 

u(x) = — k tanh(eofcc) 

f(x) = cosh~ q (kx/q) (3.22) 

where q = 1/eo- Also, from (3.21) and (3.22), we have the following exact solutions for 
and V{4>) : 

= arctan[sinh(/cx/g)] 
V((f)) = {cos0} 2(1 " 9) 
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= /c 2 sech 2 ( k / e(x) dx + c 



where the parameters q and k are related to the energy scales A and m by 

A 2 = 2g(l-g)(l/87rG) 
m 4 = A; 2 (i + 2)(l/87rG) 

Thus, f(x), 4>(x) and V(4>) solve the Einstein equations as well as the field equation for 
the real scalar <fi. Finally, the metric for the accelerated spherically symmetric domain wall 
spacetime is given by 

g = f( x ) {-At <g> dt + dx <g> dx} + (y cosh(kt)) 2 f 2 (x) {dO ® d6 + sin 2 9dip <g> d^j 

k 1 > 

(3.23) 

It is interesting to note that for a fixed x = xo, the spherical domain represented by (a 
3-dimensional De-Sitter slice) 

g\ x =x = f 2 {xo) |-dt ® dt + (i cosh(A;t)) 2 (d6' ® d6» + sin 2 6»d^ ® cty>) J 

(3.24) 

first contracts and then expands. Furthermore, if we choose a chart {r, i], y, z} with the 
following (implicit) coordinate transformations 

11 1 
r = (— ) ln(— sinh kt + — cosh kt cos #) 

K K r£ 

T] = X 

y = (cosh fct sin 8 cos ■?/>) / (sinh kt + cosh /rf cos 9) 

z = (cosh kt sin 6* sin ^) / (sinh kt + cosh kt cos #) (3.25) 
the wall spacetime (3.23) becomes (locally) plane symmetric: 

g = f 2 (i]) {-dr <g> dr + dr? <g> d?]} + exp(2kr) f 2 (i]) {dy <g> dy + dz dz} 

(3.26) 

The metric (3.23) or (3.26) is only one of many possible domain wall solutions and the 
equations (3.19)-(3.20) show that admissible choices for e(x) generate additional solutions. 
For example, the smooth function e(x) = 2 [1 — kx exp(— k 2 x 2 )] > 1 has aymptotic limits, 
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lim^i^oo e(x) = 2. Thus, e(x) belongs to the class of bounded smooth functions so that (3.21) 
is satisfied, and hence, e(x) must lead to a domain wall solution to the Einstein equations. 
The equation (3.19) for u(x) then becomes 

u'(x) 



k 2 — u{x) 2 



= -2 



- kxexp (-k 2 x 2 )] . (3.27) 

By integrating (3.27), we get 

u(x) = -k tanh [2kx + exp(-A; V) - C] (3.28) 

where C is a constant of integration which we will set to zero. The constant C determines 
where u(x) crosses the x-axis. By adjusting C the function u(x) can be made to cross the 
x-axis at the origin. Since u = f'/f we can now solve for f(x) which gives 



f(x) = exp 



-k J tanh (2kx + exp(— k 2 x 2 )^j dx 



(3.29) 



The integral in (3.29) which can not be evaluated analytically, is easily done numerically. A 
plot of f(x) is shown in Fig 1. 

The energy density p{x) and the pressure v{x) can be computed from the Einstein 
equations (3.14) and (3.15) and are shown in Fig 2 and Fig 3. The equation of state can 
no longer be found analytically but it is possible to parametricaly plot the energy density 
and the pressure as a function of position. This is shown in Fig 4. The non-zero area inclosed 
by the p - v diagram is due to the lack of reflection symmetry in the stress tensor. The kink 
in the curve is due to the bump in the potential V((f>) which is shown in figure 5. The field 
which generates V(<f>) is shown in figure 6. 

The equations (3.19)-(3.20) clearly allow one to generate two broad classes of accelerated 
spherically symmetric domain wall spacetimes but all of these solutions lead to solition-like 
configurations for the scalar field and have geodesies with the same general features. These 
solutions may or may not be reflection symmetric in x-coordinate. In the next section, we 
will show that static wall spacetimes cannot be reflection symmetric, leading to non-trivial 
asymptotic vacuum structures. 

12 



IV. STATIC PLANE SYMMETRIC DOMAIN WALLS 



We now look for static solutions to the Einstein-Scalar equations. Imposing plane symmetry, 



the metric for a static spacetime takes the following form [fLq| : 

g = f 2 (x){-dt ®dt + dx(& dx} + h 2 (x){dy <g> dy + dz <g> dz} (4.1) 
Then the Einstein's equations corresponding to the domain wall stress tensor are given by 

u' + 2uv = -(4nG)f(p- v), (4.2) 

u + 2(v' + v 2 - uv) = -(AnG)f 2 (3p + u) (4.3) 

v' + 2v 2 = -{47cG)f{p-v) (4.4) 
where, we set u = (/'//) an d v = (h'/h). Equations (4.2)-(4.4) also imply [see Appendix] 

(u' + 2v') < (4.5) 

v' = -(p + v)(u + 2v) (4.6) 

The above equations lead to a number of general properties of the plane symmetric static 
domain wall, assuming that (p + v) > 0, limi^i^oo v = and lirnui^oo p = 0. 

4.1 Negative 'Pressure' (is) : 

From the equation (4.5), it follows that (u+2v) has at most one zero. Also, for a nowhere 
zero smooth function, u, there exists an x such that v'ix) = (See Appendix, Lemma 2, 
for a proof of this statement). Then, (4.6) implies that each of the functions (u + 2v) and 
v 1 has a unique zero at x — x. Furthermore, differentiating (4.6) and using (4.5), we have 
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v"{x) = — (p + u)[u'(x) + 2v'(x)] > - which means that the unique extremum of v over R 
is a minimum at x — x. However, v satisfies the asymptotic condition lim^i^oo v — 0, and 
hence v must be negative. 

4.2 Plane Symmetric Static Wall is Not Reflection Symmetric : 

If the metric (4.1) is reflection symmetric about x = (say) then f(x) = f(—x) and hence, 
f'(x) = —f'(—x) and /'(0) = 0. Similar results hold for h(x), and u(0) = = v(0). Also, 
from the Einstein's equations (4.2) and (4.3), we have 

(y'-u') + 2v{v-u) = 

Integrating once gives, 

[v(x) -u(x)]h 2 (x) = K/2 (4.7) 

for all x where, K is a constant. For K ^ 0, equation (4.7) contradicts reflection symmetry 
of f(x) and h(x), and hence a plane symmetric static wall can not have reflection symmetry. 
If K — 0, is proportional to f(x) and with the linear changes in coordinates we have 
f(x) = h(x) and u(x) = v(x). Then from Einstein's equations, 

u' + 2u 2 = -{4nG)f{p-v) 

u + 2(v' + v 2 - uv) = -(47rG)/ 2 (3p + v) 

we find 

< 3m 2 = (8nG)f 2 u < 

since, v < for non-zero v. Thus, for K = we have v — 0, and hence p = producing a 
vacuum. This conclusion also means that a conformally flat plane symmetric static metric 
leads to a vacuum spacetime under the asymptotic conditions on v and p. 
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4.3 Asymptotic Limits of The Einstein Equations : 

If we assume, again, that pressure (u) and energy density (p) vanish in the region far from 
a thick membrane, then the following equations 

lim v = ; lim p = 

\x\— +oo \ x \— > °° 

generate the vacuum Einstein's equations: 

u +2uv = Q (4.8) 

u' + 2(y' + v 2 -uv) = (4.9) 

u' + 2u 2 = (4.10) 
From these equations, we find an algebraic relation 

(2m + v)v = 

(a) If (h'/h) = v — 0, then (/'//)' = «' = 0. In this case, we may take = 1, and 
on integration f(x) = exp(px), where p is a constant. Inserting the following coordinate 
transformations 

t — - exp(pt) sinh(jxc); x = - exp(pt) cosh(pa;); y = y; z = z 

in (4.1), we have the Minkowski spacetime: 

QMink = —dt <8> dt + dx <S> dx + dy ® dy + dz <g> dz 

(b) If 2(f'/f) + (h'/h) = 2-u + w = 0, by integrating once we have f 2 (x)h(x) = C , where C 
is a constant. Also, (4.10) implies /i 2 (x) = x, and choosing C = 1 we get / 2 (x) = \jyfx. In 
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this case, we have Taub's plane symmetric vaccum spacetime [13]: 



dTaub = —^{—dt (g> dt + dx (g> dx} + x{dy ® dy + dz ® dz\ 



4.4 Solutions to Einstein's Equations : 

In order to obtain a plane symmetric static wall solution to Einstein's equations (4.2), (4.3) 
and (4.4) we now introduce a coordinate transformation suggested by the equation (4.7) : 

1 dhix) 1 df(x) K , . 

+ (4-11) 



h(x) dx f(x) dx 2h?{x) 
If we define the non-singular coordinate change £ = £(x) by 



h 2 (x 

then considering f(x) and h(x) as functions of £, we have 



^=#T ( 4 - 12 ) 



ldh ldf K lA . 
= — - H (4.13) 

hdi fdC 2 1 1 

Since f(x)dx = f(x(£,))h 2 (x(£,))d£,, we define 

f(x(0) = F(0 and iV(0 ee /(x(0)/> 2 (x(0) (4.14) 

Thus, h 2 (x(^)) = and hence, in the new coordinate system {t, £, y } z} } the plane 

symmetric static metric (4.1) takes the following form: 

g = -F 2 (£)dt ®dt + N 2 {£)d£ ®d£ + ^$-{dy <8> dy + dz ® dz} (4.15) 

^ I? J 

Denoting differentiation with respect to £ by a ('), the equation (4.13) implies 

Einstein's equations corresponding to the metric (4.15) are given by 
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-{^G)N 2 {p- v) 



(4.17) 



3 (F^ 
2 \ F 



'F'" 



'.AT 
N 



1 /rr 

2 I N 



-(47rG)iV 2 (3p + z/) (4.18) 



1 

2 



o Kf - hr =-(4vrG)iV 2 (p-z,) 



(4.19) 



From (4.17)-(4.19) and (4.16), we have 



'F'V 



F' 



F' 



2 - + 3 r-\ + 2r-\K + K'/A\ = (8nG)N 2 p(0 



(8nG)N 2 u(0 



(4.20) 



If we set r(£) = F'(£)/F(£) + K/3 and a = -fC/6, then the Einstein equations (4.20) becomes 



2r - 3(r 2 - a 2 ) = -{8nGN 2 )p 



-3(r 2 - a 2 ) = -{8nGN 2 )u 



(4.21) 



From this form of the Einstein equations, we have 



r'(0 



a 2 - r(£Y 



(4.22) 



where the nowhere zero smooth function er(£) is given by 



3p(0 

KOI 



+ 1 



> 3 



(4.23) 



since p(£) > |^(£)l- Since v < 0, we have an equation of state in the form = — ?(£) p(£)> 



where the relation | 



3p(Q 

KOI 



+ 1 



cr(£) > 3 is equivalent to < q(£) < 1. Following 
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the procedure developed in subsection 3.2, for generating solutions to the Einstein-scalar 
equations, we obtain the simplest static domain wall solution to (4.20)-(4.21) when q is a 
constant : 



L 

In F = In 

3 



cosh f — - 
I 2L 



^ (4.24) 



where L is related to q by q = L/{2 — L) and hence, < L < 1. Now, using the equation 
(4.16) we find the metric functions: 

F 2 = (cosh(i^/2L)r^ exp ( 4 - 25 ) 
N 2 = (cosh(fT£/2L))~ 2L (4.26) 

Q = (cosh(i^/2L))-^ exp (^J (4.27) 

Equations (4.25)-(4.27) represent a two-parameter (X, L) family of plane symmetric static 
spacetimes dominated by domain walls. The 'density' (p) and 'pressure' (v) functions are 
obtained from the Einstein's equations [4.20-4.21]: 

A' 2 (2-L) sech(ff) 2 " 21 
"«» 96 G J (428) 

"® = 96^" < 429) 

The explicit forms of the (dimensionless) scalar field and the potential are given by 

0(0 = arctan (sinh J J ( 4 -30) 



V{<j)) = (cos0) 2 " 2L 
18 



(4.31) 



The energy scales, A and m, are related to the parameters, K and L, by 

mi = ^r, < 433 > 

It should be noted from (4.28)-(4.29) that the matter density and pressure are reflection 
symmetric, while the spacetime [determined by (4.25)-(4.27)] is not. The absence of reflec- 
tion symmetry of the static wall spacetime is also demonstrated by the existence of different 
asymptotic vacua defined by 

lim v — ; lim p = 

For £ — > oo, 

F 2 -> exp (-iff) (4.34) 

N 2 -> exp (-K£) (4.35) 

I — 1 — > constant (4.36) 

Comparing our results in subection 4.3, it is clear that the limits (4.34)-(4.36) implies a 
Minkowski vacuum. For £ — > — oo, 

F 2 - exp (-j^) (4.37) 
N 2 -> exp (iff) (4.38) 
?L) _> exp f (4-39) 



F7 " V3 

The fact that the limits (4.37)-(4.39) yields the Taub vacuum, can be seen as follow: 

1 2 

5 ,^^_ 00 = - exp ( — K£)dt ®dt + exp {K^)d^ <g> d£ + exp (-K£){dy ® dy + dz <g> dz} 
3 3 

(4.40) 
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Then, the coordinate transformations given by 

* = (2K/3)h 

e = (3/2A-)ln{(2tf/3)$&} 
y = (3/2K)^y 

z = (3/2K)h a 



(4.41) 
(4.42) 
(4.43) 
(4.44) 



lead to the metric for the Taub vacuum: 



g Taub = —={-dt Q ® dt + d^ <g) d£ } + t; {dy <g) dy Q + dz a ® dz Q } (4.45) 



4.5 More Static Solutions 



In (4.22), if we take a nowhere zero bounded smooth function 



<7(£) = 91 + exp(-36a 2 £ 2 ) 



which satisfies (4.23), then we have 



r(£) = — atanh 



(9a)£+(-)v^Erf(6aO-<? 



(4.46) 



A plot of r(^) is shown, for C = and a = 1/6 [or equivalently, K — 1], in figure 7. The 
function can be found by numerical integration. A plot of F(£) is shown in figure 8. 
Equation 4.16 can be integrated to give iV(£). A plot of N(£) is shown in figure 9. 



V. PARTICLE MOTION 



For domain wall spacetimes gravity is, on average, repulsive. This can be seen from 
the following argument. The condition, R\c(v,v) > for all observers v, expresses the 
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emperical fact that, on average, gravity attracts [pl|. This condition is equivalent to the 
strong energy condition, T(v,v) > — ^ for all observers v, where T = trace(T). Thus, 
violation of the strong energy condition for the Einstein-Scalar stress tensor implies that, 
on average, gravity is repulsive in the domain wall spacetime. 



A. Accelerated Spacetimes 



To obtain the geodesic eqations for a particle (of mass m) in the accelerated domain wall 
spacetime [equation (3.23)] , we consider its worldline in the equatorial plane [6 = ir/2] : 

7 : r i-f (t = t(r), x = (r), 9 = tt/2, ip = ip(r)) (5.1) 

If 7 is a geodesic (and r is an affine parameter along 7), then the 4- velocity of the particle 
is given by 

*-(#MIM3)* ™ 

with the normalisation #(7*, 7*) = — m 2 : 

~ / 2 (f) + f (£)' + f^tlkf = -m- (5.3) 

Then the geodesic 7 is specified by 

(S)'-(£)GM (5 - 4) 



cfy _ / J_\ / a 3 k 2 ' 
^ ~ l/V Uosh 2 fci 



(5.6) 
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where ao and 03 are constants of integration. For a massive particle (m = 1), 3- velocity with 
respect to an observer Z = (1/ f)d t is given by 



For any given domain wall solution these equations can be solved numerically. We will 
work with f(x) = cosh -9 (kx / q) [(3.22)] solution but the essential features of particle motion 
remain unchanged in other spacetimes obtained from (3.19)-(3.21). Since f(x) decreases 
monotonically away from x = (where it attains its unique maximum), the 3- velocity of 
the partcle with respect to the observer Z, will increase away from x = 0. This shows the 
repulsive nature of the domain wall gravitational field 

Since domain wall spacetimes allow repulsive gravity it is clear that "turning point" 
solutions in which a incoming particle is repelled by the domain wall are allowed. The phase 
space plot for one such geodesic is shown in figure 10. Since the spacetime is expanding and 
particles are accelerated away from the wall the coordinate time t along a geodesic is red 
shifted. A plot of t(r) is shown in figure 11. 



To analyse particle motion in a static domain wall spacetime we now consider the geodesic 
equations for a test particle (of mass m = 0, 1) . Such a particle may be represented by a 
curve on (M 4 , g) [with g given by (4.15) and (4.25)-(4.27)] : 




(5.7) 



B. Static Spacetimes 



7 : r 1— > 



(t = t(T)^ = ^r),y = y(T),z = z(r)) 



(5.8) 
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If 7 is a geodesic (and r is an affine parameter along 7), then the 4- velocity of the particle 



is 



with the normalisation #(7*, 7*) = — m 2 : 



-*(s)'+*(#)'+w*{@Hs)">-^ (510) 

The geodesic equations are given by 

s K ) - 5 {-^' (lh {N2y (£)' + WFy &J + WFY @0] <» 2 > 

%L = C 2 /(N/F) (5.13) 
- = C 3 /(N/F) (5.14) 



where C , C 2 , and C 3 are constant and 'prime' represents differentiation with respect to £. 

For a massive particle (m = 1), the 3- velocity with respect to an observer Z = (1/ F)d t is 
given by 

v 2 = l-(F/C ) 2 (5.15) 

Since F is a monotonically decreasing positive function, a massive particle that starts on 
the left side with a negative velocity will move away from the wall, reach a turning point 
and be accelerated back towards the wall. After the particle has passed through the wall 
it will be repelled from the wall. Thus, on the left side of the domain wall (Taub- vacuum), 
gravity is attractive while, on the right side (Minkowsky- vacuum), gravity is repulsive. A 
phase space plot for such a particle is shown in figure 12. 
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For a photon traveling along £-axis, the corresponding null geodesic is represented by 

7 : A h-> (i = t(A), f = £(A), y = a, z = 6) (5.16) 

where, A ia an affine parameter, and a, b are constants. This null geodesic is characterised 
by 

£ = 2 2 (5.18) 

9(7 „ 7 .).0^-^g) 2 + ^(|) 2 (5,9) 

Now, with respect to the observer field Z = (1/F)d t , we compute the frequency ratio for 
the photons traveling along £-axis: 

j -f = {-g(^, Z e )}/{-g(^, Z Q )} = FfoyFfo) (5.20) 

Jo 

where e and o are emission-point and observation-point respectively. From the above equa- 
tion, we note that photons approaching the wall from the left-hand (Taub) side will be 
continuously blue-shifted as they pass through the wall toward the Minkowski vacuum. 
Photons that approach the wall from the right-hand (Minkowski) side of the wall will first 
be red-shifted and then blue-shifted when they pass through the wall into the Taub side. 

The geodesic equations with respect to the generalised static wall metric discussed in 
section 4.5, can be obtained by replacing F in the above equations. The essential charac- 
teristics of particle motion are the same as discussed above. 
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VI. CONCLUSIONS 



We have found new classes of static and non-static domain wall solutions to the coupled 
Einstein-scalar field equations which are (locally) plane symmetric. The static walls are 
particularly interesting since static domains walls cannot be reflection symmetric. Conse- 
quently, static domain walls must possess different asymptotic vacua. In principle, it should 
be possible to detect static domain walls by studing the red-shift of photons. The fact that 
photons approaching the wall from the Taub side are continuously blue-shifted as they pass 
through the wall toward the Minkowski vacuum, while photons that approach the wall from 
the Minkowski side of the wall are first red-shifted and then blue-shifted when they pass 
through the wall into the Taub side provides a unique experimental signature. 

It is generally believed that small perturbations in the mass distribution in the early 



Universe could subsequently grow [20,21 1 to form galaxies and the structures we currently 



see around us. However, one of the major problems in modern cosmology is to identify 



the origin of these perturbations. The domain walls discussed in this paper could be possible 
sources of the density fluctuations through their specific gravitational interactions with the 
ambient matter. 
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APPENDIX 
Lemma 1: (u' + 2v r ) < [Equation (4.5)]. 
Proof: Equating (4.2) and (4.4) leads to 

2{v 2 — uv) — u' — v' 

Inserting the above relation in (4.3) gives 

u' + v '/2 = -(27rG)/ 2 (3p + v) < 

This equation, together with v' < [which follows from (4.4)], proves the lemma. 

Lemma 2: Let v(x) be a smooth nowhere-zero function on R. If lim^i^oo v — 0, then there 
is an x such that u'(x) = 0. 

Proof: Suppose that v' ^ for all x E M. Since v(x) is smooth, v'(x) is continuous and 
hence, either v'ix) > or v'ix) < 0. If v'ix) > 0, then v is a strictly increasing function 
on M, and hence it can not satisfy the boundary condition Xm^^^^v = 0. For v'{x) < 0, 
a similar argument leads to the same contradiction. Thus, there exists aniGl such that 
v\x) = 0. 
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FIGURES 



f( x) 




FIG. 1. Plot of f(x) as defined by Eq. 3.29 versus x. Note that the domain wall is no lonj 
centered at x = 0. 

P(x) 




FIG. 2. Plot of the energy density p{x) versus x for k = 1. 



V(x) 




FIG. 3. Plot of the pressure v{x) versus x for k = 1. Note that the pressure is negative, as 
must be for a domain wall solution. 
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FIG. 4. Plot of the energy density p(x) versus the pressure v{x) for k = 1 . 

V(4> (x) ) 
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FIG. 5. Plot of the potential V(4>) versus x for k = 1. 
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FIG. 6. Plot of the field 4>{x) versus x for k = 1. 
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FIG. 7. Plot of r(£) as defined by Eq. 4.46 versus x for K = 1. 
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FIG. 10. Phase space plot of a turning point solution for a particle with initial conditions 
x(0) = 20, x (0) = -0.5, y (0) = 0.5, y(0) = 0, z (0) = 0, z(0) = 0, t (0) = 1, t(0) = and f given by 
3.22. 




2 4 6 8 

FIG. 11. A plot of r(t) showing the redshift of the t for a particle with initial conditions 
s(0) = 1/2,0/(0) = -0.5,y'(0) = 0,y(0) = 0,z'(0) = 0,z(0) = 0,t'(0) = l,t(0) = and f given by 
3.22. 
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FIG. 12. Phase space plot of a turning point solution for a unit mass particle with initial 
conditions x(0) = 2,x'(0) = -0.5,y'(0) = 0.9188, y(0) = 0,z'(0) = 0,z(0) = 0,t'(0) = l,t(0) = 0. 
On the left hand side of the domain wall gravity is attractive, while on the right hand side of the 
domain wall gravity is repulsive. The space-time has parameters m = 1,L = 1/2, and k = 1/2 and 
metric given by g in (4.15 and 4.25-4.27) 
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FIG. 13. The red shift of r(t) along a gedodesic for a unit mass particle with initial conditions 
x(0) = 2,x'(0) = -0.5,y'(0) = 0.9188, y(0) = 0,z'(0) = 0,z(0) = 0,t'(0) = l,t(0) = 0. The 
space-time has parameters m = 1,L = 1/2, and k = 1/2 and metric given by g in (4.15 and 
4.25-4.27 ) 
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